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1 Introduction 

An important motivation for the development of the quantum cosmology programme has been to 
explain the initial conditions for the emergence of the Universe as a classical outcome. In principle 
one must find the form of the wave function ^> satisfying the Wheeler-DeWitt equation H. This 
equation describes the annihilation of the wave function by the Hamiltonian operator and since it 
admits an infinite number of solutions, one must also choose the boundary conditions in order to 
specify the wave function uniquely. Such boundary conditions must be viewed as an additional 
physical law since, by definition, there is nothing external to the Universe. In practice one assumes, 
at least implicitly, that a finite subset of all possible boundary conditions is favoured by cosmological 
observations, in the sense that the wave functions corresponding to such boundary conditions predict 
outcomes which are compatible with observations. For example, if one believes in the inflationary 
scenario, the requirement that sufficient inflation occurred, in order to solve the assorted problems of 
the standard big bang model can, in principle, restrict the number of plausible boundary conditions. 

Among the set of all possible choices the Vilenkin, or tunneling from nothing, boundary condi- 
tion ^ ||] and the Hartle-Hawking, or no-boundary, boundary condition Q have been the subject of 
intense discussion. Given the non-uniqueness of such conditions, the question arises as to the con- 
sequences of choosing different boundary conditions for the resulting wave function of the Universe 
and its corresponding probability measures. An important study in this regard is due to Vilenkin jj| , 
who considered the effects of the above boundary conditions within the context of Einstein gravity 
minimally coupled to a self-interacting scalar field. He restricted his analysis to the minisuperspace 
corresponding to the spatially closed, isotropic and homogeneous Friedmann-Lemaitre-Robertson- 
Walker (FLRW) Universe and showed that the tunneling wave function predicts initial states that 
are likely to lead to sufficient inflation, whereas the Hartle-Hawking wave function does not. 

It is sometimes argued that this result indicates that observations favour the tunneling as op- 
posed to the no-boundary boundary condition. However, the precise relation between the boundary 
conditions and the observations is determined by the specific models employed and since such mod- 
els always involve idealisations in the form of a set of simplifying assumptions, it follows that the 
above conclusion can not be made a priori. Indeed it only makes sense in general if the correspon- 
dence between the observations and the boundary conditions is robust under physically motivated 
perturbations to the underlying quantum cosmological model. 

Consequently, it is important to consider the 'stability' of the above conclusions. In particular, are 
the conclusions robust under higher-order perturbations to the Einstein-Hilbert action? Quadratic 
and higher-order terms in the Ricmann curvature tensor and its traces appear in the low-energy limit 
of superstrings ^ and they also arise when the usual perturbation expansion is applied to General 
Relativity || f?J. Such terms diverge as the initial singularity is approached, but can in principle be 
eliminated if higher-order corrections are included in the action. In four-dimensional space-times the 
Hirzebrucht signature and Euler number imply that the most general, four-dimensional gravitational 
action to quadratic order is 

S = J ' otx S~g ^2 [ -lC al3lS C af} ^ + e 1 R 2 ] , (1) 

where R is the Ricci curvature scalar of the space-time with metric tensor g^ v , g = detg^, C a( g 7 <5 
is the Weyl conformal curvature tensor, k 2 is the gravitational coupling constant and e\ and 7 are 
coupling constants of dimension (length) 2 . The action simplifies further for spatially homogeneous 
and isotropic four-geometries, since the conformal flatness of these space-times implies that the Weyl 
tensor vanishes. The effects of including quadratic terms have been investigated in Refs. |||], [To| - 
In particular Mijic et al [[l0| studied the effects of such perturbations on Vilenkin's result |J and 
found that those results remain robust in the sense that the inflationary scenario still favours the 
tunneling boundary condition in the presence of quadratic terms in the action. On the other hand 
Biswas and Guha have recently arrived at the opposite conclusion B. 
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The renormalisation of higher loop contributions introduces terms into the effective action that 
are higher than quadratic order. Consequently it is important to also study the effects of these 
additional terms. In this paper we shall investigate what happens to the wave function if an R 3 - 
contribution is present. By employing the conformal equivalence of higher-order gravity theories 
with Einstein gravity coupled to matter fields, we argue that this term represents a more general 
perturbation to the Einstein-Hilbert action than the i? 2 -correction, at least within the context of 
four-dimensional FLRW space-times. We then consider the conditional probability that an infla- 
tionary epoch of sufficient duration can occur. We estimate how the qualitative behaviour of this 
quantity changes when higher-order perturbations to the action are included. Our main result is that 
for the i? 3 -theory there exists a finite region of parameter space in which neither of the boundary 
conditions discussed above predict an epoch of inflationary expansion that leads to the observed 
Universe. We use (dimensionless) Planckian units defined hy Ti = c = G= \ throughout and define 
k 2 = 8tt. 



2 Higher-Order Lagrangians as Einstein Gravity plus Matter 

The wave function of the Universe in higher-order Lagrangian theories can be determined in one of 
two ways. It is well known that theories with a Lagrangian given by a differentiable function of the 
Ricci curvature scalar are conformally equivalent to Einstein gravity with a matter sector containing 
a minimally coupled, self- interacting scalar field [T^j , The precise form of the self-interaction 
is uniquely determined by the higher-derivative metric terms in the field equations. It follows that 
one can start either from the original action or the conformal action and derive the corresponding 
Wheeler-DeWitt equation . One takes the related Lagrangian as the defining feature of the theory 
and then applies the canonical quantisation rules. The advantage of the conformal transformation is 
that it allows the known results from Einstein gravity to be carried over to the higher-order examples 
and we shall follow such an approach in this paper. 

Consider the general, D-dimensional, vacuum theory 



d u x [ f(R) 



(2) 



where the Lagrangian f(R) is some arbitrary differentiable function of the Ricci curvature scalar 
satisfying {f(R), df (R) / dR} > and is the determinant of the D-dimensional space-time metric 
9d^v If we perform the conformal transformation fl^| 



2/(D-2) 



(3) 



and define a new scalar field 



D- 1 
D-2 



1/2 



In 



2k' 



( df(R) 
\ dR 



(4) 



the conformally transformed action takes the Einstein-Hilbert form 



S = 



fJD 



(5) 



where the self-interaction potential is given by 



U(4>) = 2k 



df[R($)\ 
dR 



-D/{D-2) 



(6) 
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Definition (||) yields a correspondence between the values of the Ricci curvature scalar R and the 
values of the scalar field <fi. We shall consider the quadratic and cubic Lagrangians 



MR) = ^(i? + eii? 2 ) 

/s(fl) = ^ 1 (R + e 1 R 2 + e 2 R 3 ) 



(7) 
(8) 



in four dimensions, where the parameters t\ and e 2 have dimensions (length) 2 and (length) 4 respec- 
tively before the introduction of Planckian units. The corresponding potentials for positive t\ and 
e 2 are given by JlTJ, : 



UfM 



1 



in 2 ei 



27k 2 e 2 2 



1 — exp(— \/2/3 k q 
exp(-2 ^/2/3k4>) 



-1 



3£2 

ei 2 



1 — exp(-y/2/3K</)) 



9e 2 
2 £l 2 
3/2 ' 



1 — exp(-\/2/3 k 



(9) 



(10) 



and are semi-positive definite for all values of c 



Figures la &: lb 



In the classical i? 3 -theory the requirement that the inflationary epoch lasts sufficiently long 
implies that the coupling constants must satisfy |e 2 | <C ei 2 0. Moreover, the observed isotropy of 
the cosmic microwave background radiation requires that ei « 10 11 JToj] . In view of these constraints 
we specify e± = 10 11 in the subsequent numerical calculations. Figures la and lb illustrate the 
behaviour of the potentials (0) and (|icj) for e\ w 10 11 and e 2 ~ 10 20 . The effect of decreasing 
the value of the parameter t\ is to increase the height of the plateau and the relative maximum of 
the potentials in the quadratic and cubic cases respectively. This reflects the fact that decreasing 
this parameter is equivalent to increasing the energy scales involved. In this sense there exists 
no continuous transformation from an i? 2 -theory to the ordinary Einstein-Hilbert action as this 
parameter approaches zero. In the neighbourhood of the origin of <j) corresponding to smaller values 
of R the quadratic term in the action dominates and the potentials in this region are equivalent. 
This can be seen by expanding the last of the three terms in the square brackets of Eq. (|h]) . The 
first-order contribution cancels the remaining terms in U f 3 and the second-order term reduces the 
form of U f 3 to that of Uf 2 . Hence the two potentials are effectively identical if the third- and higher- 
order terms in the expansion can be neglected. It is straightforward to show that this is a consistent 
approximation if 



K (f> < K 4>Umit = 




For polynomial Lagrangians with f(R) = (Y^k=i e k-iR k ) /2k 2 , the detailed form of the cor- 
responding potential U{<f) is extremely complicated and generally not expressible in an analytically 
closed form. Nevertheless, one can determine the qualitative behaviour of the potential at small and 
large (f>. Close to the origin the quadratic term in the action again dominates and the potential in 
this region is therefore similar to Eq. (^). The asymptotic behaviour at infinity, however, depends 
critically upon the combination of the highest degree n of the polynomial and the dimensionality 
D of the space-time [Q. More precisely, for D > 2n the potential is unbounded from above, for 
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D = 2n it flatens into a plateau and for D < 2n the potential has an exponentially decaying tail 
|L4j . In particular, if D < 2n the effective scalar field potential U(4>) is qualitatively equivalent to 
the cubic potential ([To|). As a result, when D — 4 the qualitative behaviour of U((j>) does not change 
relative to the cubic case as terms with n > 3 are considered, although the relative position of the 
maximum of U((j>) will be n-dependent. This implies that the n = 2 contribution is rather special in 
four dimensions, whereas the i? 3 -term is in fact a more generic perturbation. Thus, it is instructive 
to consider this case further. 



3 Behaviour of the Wave Function 



Within the context of the spatially closed FLRW minisuperspace, the Wheeler-DeWitt equation 
derived from theory (|J) has been solved for an arbitrary potential, subject to the condition that the 
momentum operator for the scalar field can be neglected (|, [|. This is self-consistent if \dV/d(j>\ <C 
max{|V|, a -2 }, where a represents the cosmological scale factor and 
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(12) 



The WKB approximations of the wave functions satisfying the quantum tunneling boundary 
condition (\Ev) and the Hartle-Hawking no-boundary proposal (^hh) then take the forms p] 



IV = (1 - a 2 V)- 1/4 exp 



HH 



= (l-a 2 V) 



-1/4 



exp 



(1 - a 2 Vf' 2 - 1 

w 

1- (1-a 2 Vf' 2 
W 



in the classically forbidden (Euclidian signature) region defined by a 2 V < 1, and 



= e l7!/i (a 2 V-l)- 1/4 
^hh = 2{a 2 V -iy 1/4 exp 





1 " 




exp 




exp 



" 1 " 




3V 


cos 





(a 2 V - If/ 2 

4 3^ 

{a 2 V - If/ 2 _ 7T 
3^ _ 4 



(13) 
(14) 

(15) 
(16) 



in the classically allowed (Lorentzian signature) region a 2 V > 1. Substituting for V(4>) from the 
potentials of the quadratic and cubic Lagrangians of Section q], it can readily be seen that the wave 
functions corresponding to the quadratic and cubic theories have very different types of behaviour, 
at least for large 4>. In the quadratic case both and \& hh remain bounded. However, for the 
cubic case ^ hh becomes divergent in the classically allowed region whilst remains regular. In 
this sense then the qualitative behaviour of the wave function satisfying the no-boundary proposal 
is fragile with respect to cubic perturbations to the action. This is significant because often the 
quadratic corrections to the action are taken as representative of higher-order perturbations. 

To proceed it is important to ensure that for the regimes under consideration the conformal 
transformation (||) remains non-singular. This is the case if the condition df(R)/dR ^ is valid for 
all values of R. The conformal transformation is singular at the point 



R 



1 

2^ 



(17) 



in the i? 2 -theory and at the point 
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for the i? 3 -theory. Since t\ and €2 are taken to be positive, these conditions imply that in both 
cases the problematic values of R lie in the region R < 0. However, for a classical, spatially closed 
FLRW model, the Ricci curvature scalar is given by 

R = 6(l-q)(~) +4, (19) 
\a J a z 

where q = —a a/a 2 defines the deceleration parameter and a dot denotes differentiation with respect 
to cosmic proper time. Now if, as is generally assumed, the Universe tunnels into the Lorentzian 
region in an inflationary phase (q < 0), it follows that R will be positive-definite. Thus, the conformal 
transformation is self-consistent in these theories. 



4 Interpretation of the Wave Function 

In the previous section we saw that the wave functions corresponding to the tunneling and the Hartle- 
Hawking boundary conditions have qualitatively different modes of behaviour for the quadratic and 
cubic theories. To see what predictive effects such changes might have, we employ the notion of 
a probability density p as is usually done. For the cases of the tunneling and the Hartle-Hawking 
boundary conditions respectively, p takes the form || 



Pv(a, 4>) — C-v exp 
PHn(a, <t>) = C HH exp 



2 



(20) 
(21) 



on surfaces of constant scale factor in the classically allowed region of minisuperspace, where the 
normalisation constants Cy and Chh are given by 



Cy- 1 



dtp exp 



C 



HH 



V(<j>)>0 



V(>p)>0 



exp 



3V(cf>) 
2 



3V(cf>) 



(22) 
(23) 



Since p{<p) is usually not normalisable, the common practice is to employ the notion of a con- 
ditional probability jl3) . One argues that the initial values of the scalar field must lie in the range 
4>min < 4>i < 4>P- The lower limit (j> m in follows from the requirement that the Universe expands at 
least until the formation of large-scale structure and the upper bound follows from the condition 
that V{<j)p) ~ 1, since the minisuperspace approximation is unlikely to be valid when the potential 
energy of the matter sector exceeds the Planck density. However, in a chaotic inflationary scenario 
there is a critical value of the scalar field, <p su f, and sufficient inflation occurs if 4>i > 4>suf but not for 
4>% < <Psuf- We must therefore calculate the conditional probability that sufficient inflation occurs 
given that fa is bounded by 4> m i n and <fip- This quantity takes the form |0 



hl t P(<t>)d<t> 



> 



^suf 



< 



< 



•) = 



(24) 



and allows us to determine which of the two boundary conditions considered here "naturally" predicts 
a phase of sufficiently long inflationary expansion. Sufficient inflation is a prediction of a theory if 
Pwl, whereas it is not if P <C 1. 

For standard reheating the minimum amount of inflation that solves the horizon problem is 
determined by the condition N = ln(a//aj) « 65, where subscripts i and / denote the values of the 
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scale factor at the onset and end of inflation respectively It is then straightforward to deduce 
from the classical field equations that 



N 



G5 



G 



<t>t 



where the value of the scalar field at the end of inflation, 0/ , is computed from the relation 



1 

12 



(25) 



(26) 



<t>=<t>f 



This condition corresponds to the breakdown of the slow-roll approximation |l6|.p| Once fat is known, 
the value of 4> S uf can be determined numerically by evaluating the integral in Eq. (^5|) . 

To understand how the probability densities (^Oj) and (^l|) change in the quadratic and cubic 
cases, we shall consider them in turn. Since (20) and ( pl| ) are usually not normalisable (unless the 
range of values that (f> can take is bounded), we set the "normalisation constants" equal to one as is 
the common practice. 



4.1 The Quadratic case 

To begin with, we note that the shape of V{fa) does not qualitatively change with changes in the 
coupling constant t\. This parameter only fixes the height of the plateau and as a result leaves the 
shapes of the two probability densities unchanged. Consequently the qualitative behaviours of the 
probability densities are robust with respect to changes in t\. Figure 2a gives a plot of pv showing 
that it starts at zero when = and asymptotically approaches a constant value. On the other 
hand, as can be seen from Figure 2b, phh decreases from infinity and asymptotically approaches 
a constant value. We should emphasise here that since the probability distribution functions ( ^p| ) 
and (|2l| ) typically take values of the order exp(±10 14 ), we, for the sake of graphical representation, 
applied non-linear scalings of the kinds pv = pv 1 ^ an d Phh = In (phh 1 ^ ) respectively (where C 
is a constant) to the two probability distribution functions. Note, however, that the values of the 
argument <fi remain uneffected by this scaling. 

Contrary to the claim of Biswas and Guha Q , the two probability distribution functions reveal 
no qualitative changes as compared to the case of "chaotic" type potentials (e.g. V(fa) — in 2 fa 2 /2) 
as discussed by Vilenkin ||| and Halliwell |H| . This means that the tunneling wave function has its 
maximum nucleation probability for the Universe coming into existence somewhere on the plateau 
of the potential V((f>), whereas the Hartle-Hawking wave function peaks near the true minimum of 
the potential at <f> = 0. Translated into initial values of the Ricci curvature scalar, this means that 
the tunneling wave function prefers values of Ri near the Planck scale, whereas the no-boundary 
wave function favours a Universe of large initial size, i.e. small Ri pi. 



Figures 2a & 2b 



We now consider the conditional probability ( p4| ) . The range of values of fa is specified by the 
range of initial values Ri . In Planckian units, where Rp — L we deduce that fao = 13.0. The value 
of far is calculated from (|2^ ) to be fat — 0.38 and condition ( |25| ) is therefore satisfied for fa u f = 2.27. 
Since the conditional probability measure (Q) essentially amounts to a comparison of areas between 
the p(fa) curve and the positive faaxia in Figures 2a and 2b, it seems obvious that the tunneling 
wave function leads to sufficient inflation whereas the no-boundary wave function does not. This is 

1 Strictly speaking, conditions (E^) and (|26|) are only valid in spatially flat FLRW models, but we are considering 
spatially closed cases in this work. However, during inflation the curvature term in the Friedmann equation is redshifted 
to zero within one Hubble expansion time and the Universe effectively becomes spatially flat at an exponentially fast 
rate. For our purposes, therefore, these expressions remain valid. 
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in line with the conclusions of Vilenkin ||| and Mijic et al ]T3] and in contrast to what is claimed by 
Biswas and Guha ||. 

4.2 The Cubic Case 

We now consider the effects of adding a cubic term to the action. In general pv is peaked around 
the maximum of V{fa) at <p max and falls off to zero on both sides. In contrast phh decreases from 
infinity near <f> = to a minimum at 4> max and diverges again as <f> — > oo. In this sense the presence 
of the cubic term drastically alters the shapes of the two probability distributions. This qualitative 
behaviour is illustrated in Figures 3a and 3b for e\ — 10 11 and e 2 = 10 20 . 

Figures 3a & 3b 

Now, regarding the location of the maximum nucleation probability, the tunneling case is un- 
ambiguous since there is only a single peak in the probability distribution function. Note, however, 
that in the cubic case this wave function favours smaller values of the initial curvature i?,; (viz. fa) 
as compared to those in the quadratic case, where they are of Planckian order. On the other hand, 
the case of the Hartle-Hawking boundary condition is ambiguous because of the presence of two 
peaks in the probability distribution function, corresponding respectively to low and high values of 

Ri- 

From a practical point of view, the question arises as to whether the Vilenkin wave function 
still predicts a phase of sufficiently long inflationary expansion immediately after tunneling into the 
Lorentzian signature region. To investigate this, we confined ourselves to the region on the left of 
the maximum in the potential (|l0|), i.e. <p < fanax- Although inflation occurs on both sides of the 
turning point, there is no end to the superluminal expansion if the field rolls down the right-hand side 
and consequently there is no reasonable mechanism of reheating [14j] . On the basis of these physical 
considerations it is therefore more appropriate to identify the upper limit </>p of the integrals in Eq. 
( pi] ) with fanax rather than with the Planck limit. 

The specific value of the conditional probability depends on the magnitude of £2 and it is therefore 
necessary to determine the relevant range of values for this parameter. We noted in Section |^ that 
£2 is bounded from above by the condition £2 <C £i 2 . As 62 is decreased relative to a fixed e\, the 
location of the maximum is shifted to larger values of (j> and eventually beyond the Planck limit <pp. 
This follows since the model reduces to the i? 2 -theory for which the potential exhibits a plateau, i.e. 
the maximum is effectively located at infinity in this case. However, according to condition ( p"l| ) the 
region over which the cubic and quadratic potentials are equivalent also increases as £2 decreases. 
The question then is whether 4>u m it grows faster or slower than fa na x ■ By explicitly calculating the 
values of 4> ma x and faimit it is found that faimit exceeds (f> ma x for all parameter values 62 < 10 20 . 
This implies that the R 2 - and i? 3 -theories are equivalent for cf> < <fi max in this range. Hence the 



results in Section 4.1 for i? 2 -theory may be carried over directly to the cubic case in this region 
of the variable fa although there is the important difference that the upper bound on fa is now 
identified with fan ax and not 4>p. 

For any given £2 the end of inflation occurs at far = 0.38 as in the i? 2 -case, since the R 3 - 
contribution is negligible at very small fa Unfortunately a direct numerical integration of Eq. (|24| ) 
can not be performed, because the integrands are typically of the orders of of exp(±10 14 ). However, 
since the probability density p is a single valued, positive-definite function of fa it follows that a 
handle on the qualitative behaviour of the conditional probability can be obtained by investigating 
how the area under the p{fa) curve changes as £2 changes. The problem then reduces to determining 
how the limits of the integrals in the numerator and denominator vary as the parameters of the 
theory are altered. 

The dependences of the parameters of interest on £2 are summarised in Table [lj. We find that 
4> su f for the potential ( |To| ) settles at the same value as in the quadratic case when £2 is of order 
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10 18 or smaller. We also find that <f> m ax rapidly approaches 4> su f in the region 10 18 < £2 < 10 20 . 
This implies that the integral in the numerator of the conditional probability ( p4| ) becomes much 
smaller than the term in the denominator for 62 > 10 18 . Consequently the Vilenkin scheme does 
not predict a phase of sufficiently long inflation in this region, contrary to the results for the R 2 - 
model. We further note that for the same range of initial values of <fi, the Hartle-Hawking wave 
function shows no qualitative change from the quadratic case. Consequently, it appears that neither 
boundary condition predicts inflation for this choice of the parameters t\ and ti- This behaviour 
occurs because the presence of the cubic perturbation severely restricts the range of initial field 
values 4n for which a phase of sufficiently long inflationary expansion is likely. 

Including the full range of values of 4>i up to the Planck limit <pp would not significantly improve 
this result in the Vilenkin scheme. In the Hartle-Hawking case, however, the integral in the nu- 
merator of ( p4|) would have a large contribution from the second peak in phh ■ However, this range 
of <f)i was excluded, as discussed above, in order to avoid the problem of exiting the inflationary 
expansion. 



£2 


1Q 20 


10 18 


10 16 


10 14 


10 12 


10 10 


10 8 


10 6 


4>p 


23.6 


21.3 


19.0 


16.7 


14.4 


13.1 


13.0 


13.0 




2.65 


4.95 


7.25 


9.56 


11.9 


14.2 


16.5 


18.8 




1.59 


2.68 


3.78 


4.94 


7.06 


9.34 


11.7 


13.0 


4>suf 


1.59 


2.24 


2.27 


2.27 


2.27 


2.27 


2.27 


2.27 



Table 1: Summarising, for different values of 62, the values of the scalar field corresponding to 
Rp = 1 (0p), the limit of <j) below which the R 2 - and i? 3 -potentials are equivalent (4>ii m it), the 
location of the maximum in the potential (cj)max) and the values of the field that just lead to sufficient 
inflation (cf> su f). We specify ei = 10 11 throughout due to microwave background considerations. As 
62 increases to order of 10 18 , the magnitudes of the quantities (j> ma x and 4> S uf become comparable 
to one another and this implies that the numerator in the conditional probability approaches zero. 
This suggests that the conditional probability will become significantly smaller than unity for values 
of e 2 > 10 18 . 



Even though the conditional probability P of Eq. (|24|) cannot be estimated numerically in this 
case, nevertheless, we present a set of values of "scaled conditional probabilities" in Appendix A 
which are obtained by applying a non-linear scaling to the probability distribution functions as 



discussed in Section 4.1. These values, which may be treated as qualitative indicators of P, also 



support the conclusions given in this section. 
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5 Discussion and Conclusions 

In this paper we have investigated how the probability of realising sufficient inflation from quantum 
cosmology is altered when higher-order corrections to the Einstein-Hilbert action are introduced. 
Our results confirm that the addition of quadratic terms to the action does not reverse the conclu- 
sions of Vilenkin j^, |3| regarding the effects of boundary conditions on the likelihood of sufficient 
inflation, in contrast to some recent claims ||. On the other hand, cubic perturbations can produce 
qualitative changes to the nature of the probability distribution function p(4>). From a physical 
point of view one is confined to consider initial values of the scalar field that allow an exit from 
the inflationary expansion. As a result the important physical (as opposed to purely mathematical) 
consequences of cubic perturbations are that they restrict the measure of allowed initial field values 
4>i that lead to sufficient inflation. This is in agreement with the classical arguments (l4|. By con- 
sidering the conditional probability (|2j) (see also Appendix A) we have argued that if the coupling- 
constant 62, which determines the strength of the i? 3 -contribution to the Lagrangian, exceeds a 
critical value, neither the tunneling nor the no-boundary boundary conditions predict an epoch of 
sufficient inflation, in the sense that the conditional probability is significantly less than unity in 
both cases. 

Our results appear to exhibit some generality in four-dimensions. As discussed in Section ||, 
the qualitative shape of the self-interaction potential V(4>) remains unaltered if general polynomial 
perturbations with a highest order term e„_i R n are considered. In general this result is true when 
D < In. This immediately implies that neither of the two probability distributions py and phh for 
the n = 3 case will be qualitatively affected under n > 3 perturbations. The qualitative conclusions 



drawn for the case of cubic perturbations in Section 4.2 therefore remain robust under higher-order 
perturbations to the action, although of course the details of what happens will depend on how the 
precise location of the maximum in the potential V((f>) is related to the highest-order term. 

However, the consequences of the quadratic and the cubic perturbations (as well as those of 
general polynomial types) depend crucially on the values of the free parameters of the system, 
namely e& (k = 1, . . . , n — 1) , D, n, as well as on the initial field values (j>i. In particular, the 
dimensionality D of the space-time is crucial in deciding the maximum degree n of perturbations 
allowed (D < 2n say) above which the perturbations would be qualitatively inconsequential, i.e. the 
system would be robust. 

Finally we remark that inflation is possible, at least at the classical level, if the field is initially 
placed to the right of the maximum in Eq. ( |l0| ) and given sufficient kinetic energy to travel over the 
hill towards </> = 0. Unfortunately, our analysis can not consider this possibility since the scalar field 
momentum operator in the Whcclcr-DeWitt equation then becomes important and the solutions 
( p0| ) and ( f2l|) are no longer valid. Furthermore, if one is prepared to include the effects of the R 3 - 
contribution in the action, the cubic term RDR should also be considered. In this case the effective 
theory resembles Einstein gravity minimally coupled to two scalar fields after a suitable conformal 
transformation on the metric []l4| and in principle a similar analysis to the one presented here can 
be followed for this more general case. We shall return to some of these questions in future. 



Acknowledgements 

We thank Jonathan J Halliwell for helpful remarks. HvE is supported by a Grant from the Drapers' 
Society at QMW. JEL is supported by the Science and Engineering Research Council (SERC), 
UK, and is supported at Fermilab by the DOE and NASA under Grant No. NAGW-2381. RT is 
supported by the SERC, UK, under Grant No. H09454. 



REFERENCES 



10 



Appendix A 



As was pointed out in Section 4.1, the integrands involved in the definition of conditional prob- 
ability typically have magnitudes of order exp(±10 14 ), which makes the numerical calculation 
of the integrals not possible in practice. Now due to the nature of these numbers no linear 
scaling of the probability function p can bypass this difficulty. The question then arises as to 
whether appropriate non- linear scalings exist which keep the conditional probability P invariant. 
To see that there do not, recall that the only scalings that leave the Wheeler-DeWitt equation, 
H ^ = 0, of the D-dimensional minisuperspace models of Quantum Cosmology invariant are given 
by H = £1~ 2 H, <J> = £! 7, I'^_ff , ]/ = W~ 2 H * = 0, (O(g) is an arbitrary function of the minisu- 
perspace co-ordinates q) provided 7 and £ (a free parameter in the Wheeler-DeWitt equation) are 
given by 7 = (2 - D)/2 and £ = —{D - 2)/8(D - 1) respectively @- Effectively this amounts to 
a redefinition of the potential U(q) and the DeWitt metric of minisuperspace f a(3 (q), which occur 
in the Hamilton operator H. More importantly, under such scale transformations the conserved 
probability current density j a defined from $ remains unchanged. This freedom, however, is not of 
much use in bypassing the numerical difficulty mentioned above in order to obtain quantitative val- 
ues for P. Nevertheless, if we confine ourselves to qualitative information, we may choose non-linear 
(but monotonic) scalings of p, which, while violating the invariance properties of the model, would 
nevertheless supply us with a qualitative indicator of P. This is not dissimilar to the way non-linear 
scalings of functions are employed for the purpose of graphical representation. 

To calculate a qualitative indicator of P we define the non-linearly scaled conditional probability 
P as 

JtL P 1/C {<t>)d4> 

where C is the index of non-linear scaling. Clearly such a scaling will not change the qualitative 
behaviour of pv and the values of its argument <f>, and therefore the values of the boundaries of the 
integrals occurring in Eq. (^4|) (as listed in Table remain the same. Furthermore, such scalings 
leave P invariant in the limiting cases where P = and P = 1. 

Here we chose C = 10 14 . Table || gives the values of P as a function of €2 for the Vilenkin wave 
function in the case of the i? 3 -model, calculated for e\ = 10 11 and the boundary values of <fi given 
in Table 0. We approximated 4>min by </>/ = 0.38. For the corresponding value of P for the Vilenkin 



P( fa > (f> suf I <f> min <<pi<<p P )= 7777— — ' ( 27 ) 



model in the i? 2 -case of Section 4.1 we found P = 0.85 



£2 


10 20 


10 18 


10 16 


10 14 


10 12 


10 10 


10 8 


10 6 


p 


0.00 


0.20 


0.45 


0.59 


0.72 


0.79 


0.84 


0.85 



Table 2: Behaviour of the non-linearly scaled c ond itional probability distribution P for the Vilenkin 
wave function in the i? 3 -model of Section 4.2. We specify ei = 10 11 throughout. 

As can be seen from Table |[ the behaviour of P supports the conclusions drawn in Section 4^2 
on the basis of qualitative analysis. 
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Figure Captions 

Figure 1: (a) The effective self-interaction potential (|^) corresponding to the i? 2 -theory with e\ = 
10 11 . The scalar field and magnitude of the potential have been rescaled via Eq. ( |T2] ) to enable 
easy comparison with the results of Section 4; (b) The rescaled effective interaction potential ( |l0| ) 
corresponding to the i? 3 -theory with t\ = 10 11 and £2 = 10 20 . 
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Figure 2: (a) The Vilenkin probability distribution pv{4>) for the i? 2 -theory with a rescaling pv(4>) — 

[ exp(— 2/3V) ] 10 ; (b) The Hartle-Hawking probability distribution phh{4>) for the i? 2 -theory 

with a rescaling phh{4>) = l n [ exp(2/3F) ] . We choose these particular rescaled values of p((f>) 
in order to obtain easily intcrpretable plots from our numerical programme. 

Figure 3: (a) The Vilenkin probability distribution pv(4>) for the i? 3 -theory with the same rescaling 
as for Figure 2a; (b) The Hartle-Hawking probability distribution phh(<P) for the i? 3 -theory with 
the same rescaling as for Figure 2b. 



